Towards T\— limited magnetic resonance imaging using Rabi beats 
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Two proof-of-principle experiments towards Ti-limited magnetic resonance imaging with NV cen- 
ters in diamond are demonstrated. First, a large number of Rabi oscillations is measured and it is 
demonstrated that the hyperfine interaction due to the NV's 14 N can be extracted from the beat- 
ing oscillations. Second, the Rabi beats under V-type microwave excitation of the three hyperfine 
manifolds is studied experimentally and described theoretically. 



INTRODUCTION 

The ever increasing need for high resolution imaging 
in the life sciences, material science, and more recently 
quantum information processing, has led over the past 
decade to the development of a variety of methods that 
try to overcome the limits set by optical diffraction [TT] . 
Exploiting the non-linear behaviour of fluorescent labels, 
a spatial resolution down to few nanometers was demon- 
strated using optical techniques such as STED 6J , PALM 
[3], and STORM [T3]. A technique whose resolution is 
independent of the wavelength, is magnetic resonance 
imaging. Magnetic resonance imaging -and nano scale 
magnetic- and electric field sensing- are key technolo- 
gies in various areas of science. Recently, magnetic reso- 
nance imaging as well as scanning probe magnetic imag- 
ing on scales relevant to molecular biological processes 
was demonstrated [2 [8]. A spin label that is suitable 
for high resolution measurements are nitrogen vacancy 
centers embedded in diamond nano crystals |17j . Ow- 
ing to their optical adressability, photostability and long 
spin coherence time recently sub lOnm spatial resolution 
as well as the detection of magnetic fields close to indi- 
vidual electron spins has been demonstrated using NV 
nano sensors [U |8]. The spatial resolution of a magnetic 
resonance measurement is determined by the strength of 
the field gradient and the effective linewidth of the elec- 
tron spin transition. Pulsed imaging modes, such as the 
Hahn-Echo and CPMG sequence can be used to decrease 
the effective linewidth. The continuous driving of Rabi 
oscillations has been proposed to decrease the effective 
linewidth down to the limit given by the spin decay time 
[16) . In this paper we perform a proof-of-principle study 
of this imaging mode. The feasibility of the method is 
demonstrated by resolving the hyperfine interaction of 
the NV center due to its 14 N nuclear spin. 

In the simplest magnetic resonance imaging mode, a 
CW electron spin resonance (ESR) measurement is per- 
formed, and the effective linewidth is given by the spin 
dephasing time T 2 *, which is typically short and only of 
the order of few micro seconds in case of the NV center. 



In the Hahn-Echo sequence, 7r/2 — it — tt/2 pulses are used 
to refocus precessing spins, such that noise generated by 
slowly fluctuating spins in the environment is canceled. 
For the Hahn-Echo sequence, the effective linewidth is 
given by the dephasing time Ti of the spin system, which 
is larger than T 2 *, and can reach several hundred mi- 
croseconds in the case of the NV center 0. More ad- 
vanced pulse sequences, such as the CPMG sequence [9], 
hold promise to decrease the effective linewidth down to 
the limit set by T[ [5Ill4j. Another elegant imaging mode 
is based on the continuous driving of Rabi oscillations. 
The decay time of Rabi oscillations is given precisely by 
T-f and depends on the microwave power. When the mi- 
crowave power is high and the Rabi frequency approaches 
the transition frequency - corresponding to the limit 
of the rotating wave approximation- Tf approaches T\. 
Among the various magnetic resonance imaging modes, 
this mode may potentially achieve the longest decay time 
and highest resolution. Over the past years there has 
been a continued effort to apply this magnetic resonance 
imaging mode to the NV center [TB] . However, this poses 
great challenges. In this mode, spatial information is ex- 
tracted from beat frequencies observed in the Rabi oscil- 
lations. A high spatial resolution requires the measure- 
ment of a large number of Rabi oscillations. Here we 
demonstrate the feasibility of this imaging mode by real- 
ising the measurement of a large number (> 500) of Rabi 
oscillations, and resolve the hyperfine splitting due to the 
NV's 14 N nuclear spin from the Rabi beats. Moreover, 
we investigate experimentally and describe theoretically 
the observed Rabi beats in case of a V-type energy level 
scheme, which is the desired excitation mode in a high 
resolution experiment. 

The principle of magnetic resonance imaging using 
Rabi beats is illustrated in Fig. [I] We envision a sample, 
such as a living cell, that is tagged with nano diamonds 
that contain single NV-centers. The NV electron spins 
are inititalized and read out optically using a confocal 
microscope. An inhomogeneous microwave field - such 
as the one generated by a coplanar waveguide - is ap- 
plied to the sample. The microwaves are matched to the 
m s = — > ±1 electron spin transition and drive Rabi 
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FIG. 1. Principle of magnetic resonance imaging based on 
Rabi beats. 



oscillations between the groundstate spin sublevels. The 
Rabi frequency depends on the microwave power, thus 
encoding the position and angle within the microwave 
field. From an a priori knowledge of the microwave field 
(or from a reference measurement), the position and an- 
gle of the individual NV-centers with respect to the elec- 
trodes can be determined. To obtain the NV position in 
two spatial directions, two successive measurements are 
performed with microwave gradients applied in different 
directions using a two dimensional waveguide structure. 

In this magnetic resonance imaging mode, the resolu- 
tion is given by the decay time T-f of the Rabi oscillations, 
the Rabi frequency f2, and the strength of the microwave 
gradient. The number N of observable Rabi oscillations 
is linked to T-f and the Rabi frequency as 



~2tT 



= N, 



(1) 



For the present coplanar waveguide, the field gradient is 
determined by the width G of the gap, and the resolution 
Sx follows approximately as 



5x = — . 
N 



(2) 



The maximum achievable Rabi frequency is of the order 
of the Lamor frequency of the spin [3] . In this case, the 
decay time Tf of the Rabi signal becomes T\ and the 
resolution approaches 



Sx = 



2irG 



(3) 



where SI = uj = D sa 2.88GHz is equal to the electron spin 
transition. For NV centers, the T\ time is typically few 
ms at room temperature |12j . such that a ratio T\juj — 
10~ 6 could in principle be achievable. With G = 10/im a 
spatial resolution of O.Olnm would be reached. However, 
this poses stringent requirements on the stability of the 
applied microwave field. Both the power stability and the 
spatial stability of the microwave field must be at least 
as good as the ratio Tf/Q, and to achieve sub Angstrom 
resolution, the microwave wire should not drift by more 
than a fraction of an atomic diameter. 
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FIG. 2. ESR spectrum of the NV center showing the hyperfine 
splitting of the m a = — ¥ —1 spin level due to 14 N. The 
spectrum was recorded with an applied DC magnetic field 
of about 150G. Detuning is denoted relative to the lowest 
hyperfine transition corresponding to 2. 4524GHz. 



RESULTS AND DISCUSSION 

To demonstrate the observation of Rabi beats, we 
study a single NV-center in a fixed microwave field. Rabi 
beat measurements with a single NV center are possible 
owing to hyperfine interaction with the NV's 14 N nuclear 
spin, that results in a splitting of the groundstate spin 
manifold into several hyperfine sublevels. Each hyperfine 
transition has a slightly different Rabi frequency, which 
results in a beat signal in the measured Rabi oscillations. 
These hyperfine beats will also be seen in high resolu- 
tion Rabi beat imaging data as a modulation, and it is 
imortant to identify and assign them correctly. In here 
we are interested in the hyperfine interaction caused by 
the 14 N nucleus (nuclear spin J = 1) of the NV cen- 
ter, which splits the m s = ±1 ground state spin levels 
each into three hyperfine sublevels with an energy split- 
ting of 6 fa 2.18MHz. This energy splitting is seen in the 
CW electron spin resonance (ESR) measurement shown 
in Figj2j In here, the ESR measurement is performed on 
the m s = — > — 1 electron spin transition. Zero detuning 
corresponds to a driving microwave field of 2. 4524GHz. 
Note, that a permanent magnetic field with a component 
along the NV axis of about 150G has been applied to 
shift the m s = — > — 1 transition away from its zero 
field value (2.88GHz) owing to the Zeeman effect, and 
thus to separate it from the m s = — > +1 spin tran- 
sition (which is shifted to higher frequency of 3.3GHz). 
This separation between the two electron spin transitions 
is necessary in order to ensure that the strong microwave 
field does not simultaneously drive both spin manifolds. 

We now consider the Rabi beats expected for the given 
hyperfine splitting. For a driving microwave field that is 
detuned by a small frequency S from the resonant tran- 
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FIG. 3. Electron spin Rabi oscillations. Solid blue line: exper- 
imental data, solid magenta line: result from Fourier trans- 
form (FigHJ. Three beating cosine with frequencies corre- 
sponding to the three local maxima in the FFT are plotted. 



sition, the Rabi frequency Q is increased as compared to 
the resonant Rabi frequency Q , following |15j 



ft = \ V& + 5 2 . 



(4) 



For small detuning (8 <§; f^o), the shift of the Rabi fre- 
quency SQ becomes 



sn = 



2Q ' 



(5) 



Suppose the microwave frequency corresponds to the 
lowest hyperfine transition (as in the present experi- 
ments). In this case, the lowest hyperfine transition 
is driven with Rabi frequency Qq, while the two other 
hyperfine transitions are driven with faster Rabi fre- 
quencies corresponding to detunings 5q — 2.18MHz and 
5_i = 4.36MHz. Using Eqfl we can evaluate the ex- 
pected Rabi frequencies. The resulting signal is the sum 
of three oscillations, whose frequencies can be extracted 
by Fourier transformation. Also, the three frequencies 
result in several beat frequencies in the Rabi signal. 

To measure Rabi beats, a coplanar waveguide 
(gap=10/xs, width=10/is) was fabricated directly onto 
a type Ha diamond sample with (100) crystallographic 
orientation and less than 4ppb nitrogen concentra- 
tion (element6, electronic grade). A DCmagnetic field 
was applied with a permanent magnet. A microwave 
field resonant to the lowest of the hyperfine transi- 
tions (2. 4524GHz) was generated with a synthesizer 
(Rhode&Schwarz SMIQ) and applied to the sample 
through a fast switch (minicircuits ZASW-2-50DR). A 
single NV center was identified in the coplanar gap us- 
ing a home built confocal fluorescence microscope. Rabi 
oscillations were measured using a standard laser- and 
microwave pulse sequence [7] . Figure [3] shows the cor- 
responding Rabi oscillations. The data shows a base 
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FIG. 4. Fourier transform of the Rabi signal. The data shows 
three distinct Rabi frequencies corresponding to the individ- 
ual hyperfine levels. 



oscillation of about 42MHz (see inset). Due to the 
beating this is approximately twice the base Rabi fre- 
quency. The two beat frequencies <5f2o = 100kHz and 
<5f2+i = 400kHz correspond to the m s = and m s = +1 
hyperfine transition. Thus in total three Rabi oscil- 
lations are observed. Their individual frequencies are 
seen in the Fourier transform (Figure [2]), that shows 
n_ x = 22.2MHz, £l = 22.3MHz, and fi +x = 22.6MHz. 
Note that the FFT data is less pronounced due to the 
relatively sparse sampling of less then 10 points per pe- 
riod (see also zoom in Fig. [3]) and power drift during 
the measurement. The sparse sampling was necessary to 
keep the total measurement time (several days for Fig. 
[3]) within reasonable bounds. To complete our measur- 
ment of the hyperfine splitting from the Rabi beats, we 
use Eq| and convert the measured beats into energy level 
shifts. We find S = 2.1MHz and 5_ x = 4.2MHz, which 
is in good agreement with the ESR spectrum (Figj2|. 

The Rabi oscillations show a decay Tj* « 25/xs. This 
decay is much much smaller than the expected T-f ss 1ms 
and is limited by microwave power drift as we show be- 
low. With the present decay T*, the resolution of the 
hyperfine measurement is as follows. The resolution of 
a hyperfine measurement from the Rabi beats is deter- 
mined by the Rabi frequency and the number N of visi- 
ble Rabi oscillations. To measure a small change in the 
Rabi period, we need to measure a number of oscillations. 
Combining Eq.(| with Eq.((]), we have 



S = 



2tt 



VA y/TTf' 



(6) 



In our case we find S ~ 6MHz, which is in good quali- 
tative agreement with our observations, since our data, 
in particular the Fourier transform, Fig|4j implies that 
we can indeed resolve the 4.32MHz detuning related to 
the farther detuned hyperfine level (22.6MHz peak in the 
FFT), however, we can just barely resolve the 2.16MHz 
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FIG. 5. Power drift. The relative change of the Rabi period 
is compared to the relative change of the microwave power. 



detuning related to the lesser detuned hyperfine level 
(22.3MHz peak in the FFT). In the present measure- 
ments the total number of Rabi oscillations is about 500, 
which translates into an equivalent spatial resolution of 
about lOnm. We expect that an improvement of the 
power stability to a level better than 10 -4 should be 
achievable without a large technical effort, and that spa- 
tial drifts of the microwave field should still be negliga- 
ble in this case. Thus a spatial resolution of about lnm 
should be achievable in a realistic measurement. 

Any fluctuation or drift of the applied microwave power 
result in a small shift of the Rabi frequency over time, 
that washes out the oscillations. To proof that the ob- 
served decay is caused by power drift, we monitor the 
power transmitted through the sample and compare it 
to the Rabi period. This measurement is realized by suc- 
cessively performing short Rabi measurements in a small 
window ranging from 25.0/xs to 25.1/Lts. The result is 
shown in Figure [5] Since the Rabi period is inversely 
proportional to the frequency and the frequency is pro- 
portional to the square root of the power, we can express 
the relation between the relative change of the Rabi pe- 
riod as 



5T 



ISP 
27' 



(7) 



The data confirms that the decay time of the Rabi oscilla- 
tions is limited by drift of the microwave power, which is 
slightly better than 10^ 3 over 24h in our case. Note that 
we also observed a change of the transmitted microwave 
power when the microscope objective was moved (about 
10~ 3 / /im). This is explained by a change of the parasitic 
capacitance of the coplanar waveguide, which plays a role 
owing to the close proximity of the high N. A. microscope 
objective to the sample (working distance 200/im). 

With the fast Rabi frequency required for a high spatial 
resolution, it is no longer feasible to shift the m s = ±1 
spin levels sufficiently far apart from each other such that 



a single spin transitions is driven selectively. This is be- 
cause the separation between the m s = ±1 levels needed 
to be larger than the Rabi frequency. This would in turn 
require a large magnetic field. However at large mag- 
netic field, level mixing causes the optical ESR contrast 
to vanish, unless the magnetic field is aligned parallel to 
the NV axis [TU]. Alignment of the magnetic field par- 
allel to the NV axis is however not possible in a sample 
with many randomly distributed NV centers. Thus, in 
a high resolution measurement, we will always measure 
with a small DC bias field and hence we will simulsta- 
neously drive the m s = ±1 spin transitions. In this case, 
six hyperfine levels are present whose Rabi oscillations 
beat with each other. This situation differs from the pre- 
vious experiments. While in the previous case, the Rabi 
signal was an incoherent sum of three oscillations (during 
a measurement, the nuclear spin switched randomly and 
slowly between the different spin states), now each nu- 
clear spin manifold forms a V-type energy level scheme 
that is driven simultaneously. In the following, we eval- 
uate the corresponding Rabi beat signal. 

The NV's spin Hamiltonian with external magnetic 
field reads 



with 



H = H 

Hzfs = D (S z - 
H Z = 5 /3BS ; 



ZFS 



(8) 

(9) 
(10) 



After transforming into the energy Eigenbasis where we 
name the m s — state as |0) and the two state that 
correspond to the m s = ±1 states as |1) and |2). Now 
the semiclassical microwave field is applied 



H 



MW 



Ae'^S,. 



(11) 



here we assume, that both transitions have the same tran- 
sition matrix elements, which breaks down with rising 
E. Now do the rotating wave approximation (RWA) and 
transform into the rotating frame which leaves us with a 
Hamiltonian of the form 



OA A 
H rot = [ A A-<5 

A A + S 



(12) 



where 5 = A_E| 1 ^| 2 )/2 and A the detuning of the mi- 
crowave from |1) + S. 

To obtain an analytical solution we set A = 0. Cal- 
culation of the Eigenvalues leads to the Rabi frequency 
\J2 V 2 + S 2 which, were the two levels detuned is \[2 
times the Rabi frequency of driving a single level. We 
now take a look at the time evolution of the system with 
initial state |0). The dynamics of |0) are 



(<5 2 + 2A 2 cos [V2A 2 + <5 2 1])' 
(2 A 2 + S 2 ) 2 



(13) 
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FIG. 6. ESR spectrum of the second NV center without DC 
bias field. The two hyperfine triplets corresponding to the 
m s — —1 and m a = +1 spin levels overlap in the central 
dip. Detuning is denoted relative to the central hyperfine 
transition corresponding to 2. 8706GHz. 
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FIG. 7. Rabi oscillations under V-type microwave excitation. 
Solid blue line: experimental data. Red and green markers 
denote the extracted beat frequencies. Solid magenta line: 
result of the extracted beat frequencies. Three beating cosine 
with the corresponding frequency shifts are plotted. 



from here we see, that there are actually two oscillation 
frequencies involved, the Rabi frequency and double the 
Rabi frequency corresponding to the cos 2 term. The lat- 
ter will be the base Rabi frequency observed in experi- 
ments, and we find the identity for the Rabi beat as 



sn = 



25* 
On 



(14) 



which is different from the previous expression, Eq.(y). 

To measure Rabi beats with V-type microwave excita- 
tion, we remove the DC bias field. Moreover we pick a 
different NV center that is closer to the central conductor 
and feels a higher microwave power. Figure [6] shows the 
ESR spectrum, with a total of five dips. The central dip 
is about twice larger as compare to the other dips. In 
this experiment, the Zeeman splitting due to a residual 
magnetic field (earth magnetic field as well as magne- 
tized parts of the setup) closely matches the hyperfine 
splitting. As a consequence, the two hyperfine triplets 
corresponding to the m s = — 1 and m s = +1 spin transi- 
tion are shifted in such a way that the highest peak of the 
lower triplet overlaps with the lowest peak of the higher 
triplet, resulting in five dips with a pronounced central 
dip. Rabi oscillations are now driven with the microwave 
frequency tuned in resonance with the central dip cor- 
responding to 2. 8706GHz. This case corresponds to the 
theoretical case A — treated above, with <$ = 2.18MHz 
and S+i = 4.36MHz for the driving of the mj = and 
rrij = +1 manifold, respectively. 

The Rabi oscillations and their Fourier transform are 
shown in Figure [7] and [HJ The first observation is that 
the Rabi frequency is increased by about a factor of two 
as compared to the previous case. Two effects contribute 
here. First, with lambda type driving the Rabi frequency 
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FIG. 8. Fourier transform of the Rabi oscillations under V- 
type microwave excitation. Peaks around 42MHz and 21MHz 
correspond to the base Rabi oscillations and a modulation 
at half the frequency. The inset shows the beating Rabi fre- 
quencies around 42MHz. Labels mark the three frequencies 
as extracted from the beat analysis of the raw data. The FFT 
resolves the fast beat, while it does not resolve the slow beat 
due to sparse sampling. 



is enhanced by a factor of as compared to driving 
a single transition, and second, the Rabi frequency is 
additionally increased due to the higher microwave power 
in the close vicinity of the central conductor. The second 
observation is a fast and a slow beat with frequency of 
about 185kHz and 812kHz, respectively. These beats are 
seen most clearly in the Rabi oscillations. Note that the 
slower beat is not resolved in the FFT (Fig(8j) due to 
power drift and sparse sampling. We can now convert 
the observed beatings to energy level shifts according to 
Eq.jjl. We find S = 2.0MHz and S +1 = 4.1MHz, which 
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is in good agreement with the ESR spectrum. 



SUMMARY AND CONCLUSIONS 

We have performed two proof-of-principle experiments 
towards T\ limited magnetic resonance imaging with NV 
centers in diamond. First, we have demonstrated the 
measurement of a large number (> 500) of Rabi flops, 
and we have shown that the hyperfine interaction due to 
14 N can be resolved from such a measurement. Second, 
we have studied the Rabi beats without a large DC bias 
field, where the nuclear spin manifolds form V-type en- 
ergy level schemes. The base Rabi frequency is increased 
by v2 and the time evolution of the state population is 
modulated by an oscillation with half the base Rabi fre- 
quency. While the present experiments were performed 
in bulk, in the future, it will be an important step to 
demonstrate Rabi beat imaging with NV centers embed- 
ded in diamond nano crystals. The present experiments 
are relevant to quantum information processing in di- 
amond. The ability to drive a large number of Rabi 
flops could allow to precisely control the spin state of 
several (detuned) NV centers simultaneously. We envi- 
sion a microwave pulse with precisely adjusted length 
that performs independent unitary transformations on a 
quantum register, such as rotate an NV A to |1) while 
simultaneously rotating an NV B to the superposition 
V / 2(|0) + |l>. 
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